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ABSTRACT 

Mkn 421 was repeatedly observed with BeppoSAX in 1997-1998. The source showed a very rich phenomenol- 
ogy, with remarkable spectral variability. This is the second of two papers presenting the results of a thorough 
temporal and spectral analysis of all the data available to us, focusing in particular on the flare of April 1998, 
which was simultaneously observed also at TeV energies. The spectral analysis and correlations are presented 
in this paper, while the data reduction and timing analysis are the content of the companion paper. The spectral 
evolution during the flare has been followed over few ks intervals, allowing us to detect for the first time the peak 
of the synchrotron component shifting to higher energies during the rising phase, and then receding. This spectral 
analysis nicely confirms the delay of the flare at the higher energies, which in Paper I we quantified as a hard 
lag of a few ks. Furthermore, at the highest energies, evidence is found of variations of the inverse Compton 
component. The spectral and temporal information obtained challenge the simplest models currently adopted for 
the (synchrotron) emission and most importantly provide clues on the particle acceleration process. A scenario 
accounting for all the observational constraints is discussed, where electrons are injected at progressively higher 
energies during the development of the flare, and the achromatic decay is ascribed to the source light crossing 
time exceeding the particle cooling timescales. 

Subject headings: galaxies: active — BL Lacertae objects: general — BL Lacertae objects: individual (Mkn 
421) — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

Blazars are radio-loud AGNs characterized by strong vari- 
ability, large and variable polarization, and high luminosity. Ra- 
dio spectra smoothly join the infrared-optical-UV ones. These 
properties are successfully interpreted in terms of synchrotron 
radiation produced in relativistic jets and beamed into our di- 
rection due to plasma moving relativistic ally close to the line 
, of sight (e.g. Urry & Padovani 1995). Many blazars are also 
strong and variable sources of GeV 7-rays, and in a few ob- 
jects, the spectrum extends up to TeV energies. The hard X to 
7-ray radiation forms a separate spectral component, with the 
luminosity peak located in the MeV-TeV range. The emission 
up to X-rays is thought to be due to synchrotron radiation from 
high energy electrons in the jet, while it is likely that 7-rays 
derive from the same electrons via inverse Compton (IC) scat- 
tering of soft (IR-UV) photons -synchrotron or ambient soft 
photons (e.g. Sikora, Begelman & Rees 1994; Ghisellini & 
Madau 1996; Ghisellini et al. 1998). The contributions of these 
two mechanisms characterize the average blazar spectral energy 
distribution (SED), which typically shows two broad peaks in 
a vF v representation (e.g. Fossati et al. 1998a). The energies 
at which the peaks occur and their relative intensity provide 
a powerful diagnostic tool to investigate the properties of the 
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emitting plasma, such as electron energies and magnetic field 
(e.g. Ghisellini et al. 1998). In X-ray bright BL Lacs (HBL, 
from High-energy-peak-BL Lacs, Padovani & Giommi 1995) 
the synchrotron maximum occurs in the soft-X-ray band. 

Variability studies constitute the most effective means to con- 
strain the emission mechanisms taking place in these sources as 
well as the geometry and modality of the energy dissipation. 

The quality and amount of X-ray data on the brightest 
sources start to allow thorough temporal analysis as function of 
energy and the characterization of the spectral evolution with 
good temporal resolution. 

Mkn 421 is the brightest HBL at X-ray and UV wavelengths 
and thus it is the best available target to study in detail the prop- 
erties of the variability of the highest frequency portion of the 
synchrotron component, which traces the changes in the energy 
range of the electron distribution which is most critically af- 
fected by the details of the acceleration and cooling processes. 

This paper is the second of two, which present the uniform 
analysis of the X-ray variability and spectral properties from 
BeppoSAX observations of Mkn 421 performed in 1997 and 
1998. In Paper I (Fossati et al. 2000) we presented the data re- 
duction and the timing analysis of the data, which revealed a re- 
markably complex phenomenology. The study of the character- 
istics of the flux variability in different energy bands shows that 
significant spectral variability is accompanying the pronounced 
changes in brightness. In particular, a more detailed analysis 
of the remarkable flare observed in 1998 revealed that: i) the 
medium energy X-rays lag the soft ones, ii) the post-flare evo- 
lution is achromatic, and iii) the light curve is symmetric in the 
softest X-ray band, and it becomes increasingly asymmetric at 
higher energies, with the decay being progressively slower that 
the rise. 

The general guidelines which we followed for the data re- 
duction and filtering are described in Paper I (in particular in §2 
and §3.2). Here we will only report on details of the treatment 
of the data specific to the spectral fitting. 
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Evolution of Mkn 421: II. Spectral Analysis and Phsyical Constraints 
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FIG. 1 . — Light curves for the 1997 (top panel) and 1998 (bottom panel) cam- 
paigns, for the 4-6 keV MECS band. The binning time is 1500 s. The reference 
time TJD=10567.0 corresponds to 1997-April-29:00:00 UT, and TJD=10924.0 
corresponds to 1998-April-21:00:00 UT. 



The paper is organized as follows. In Sections §2 and §3 
we briefly summarize the basic information on BeppoSAX and 
the 1997 and 1998 observations. The results and discussion 
relative to the spectral analysis are the content of §4 and §5. In 
particular, the observed variability behavior strongly constrains 
any possible time dependent particle acceleration prescription. 
We will therefore consider these results together with those of 
the temporal analysis, discuss which constraints are provided 
to current models and present a possible scenario to interpret 
the complex spectral and temporal findings (§5.4). Finally, we 
draw our conclusions in §6. 

2. BEPPOSAX OVERVIEW 

For an exhaustive description of the Italian/Dutch BeppoSAX 
mission we refer to Boella et al. (1997) and references therein. 
The results discussed in this paper are based on the data ob- 
tained with the Low and Medium Energy Concentrator Spec- 
trometers (LECS and MECS) and the Phoswich Detector Sys- 
tem (PDS). The LECS and MECS have imaging capabilities in 
the 0. 1-10 keV and 1 .3-10 keV energy band, respectively, with 
energy resolution of 8% at 6 keV. The PDS covers the range 
1 3-300 keV. 

The present analysis is based on the SAXDAS linearized 
event files for the LECS and the MECS experiments, together 
with appropriate background event files, as produced at the 
BeppoSAX Science Data Center (REV 0.2, 1 . 1 AND 2.0). The 
PDS data reduction was performed using the XAS software 
(Chiappetti & Dal Fiume 1997) according to the procedure de- 
scribed in Chiappetti et al. (1999). 

3. OBSERVATIONS 

Mkn 421 has been observed by BeppoSAX in the springs of 
1997 and 1998. For reference, the BeppoSAX light curves for 
the 4-6 keV band are reported in Fig. 1. Most of the forth- 
coming analysis is focused on the spectral variability observed 
during the flare of 1998 April 21 st , clearly visible in the bottom 
panel. The journal of observations is given in Table 1 of Paper I. 



4. SPECTRAL ANALYSIS 

In order to unveil the spectral variability during and after the 
flare, i.e. perform a time resolved spectral analysis, we subdi- 
vided the whole dataset in sub-intervals corresponding to sin- 
gle BeppoSAX orbits (42 for the 1997 dataset, 16 on 1998 April 
21 st , 19 on 1998 April 23 rd ) or a grouping of them sufficient to 
reach the same statistics (these are reported in Table 1, together 
with the median time [UTC] of each sub-interval). 

LECS and MECS spectra have been accumulated as de- 
scribed in Paper I. LECS data have been considered only in the 
range 0.12-3 keV due to calibration problems (a spurious hard- 
ening) at higher energies (Guainazzi 1997, private comm.) 9 

The nominal full resolution spectra (i.e. channels #1 1-285 
for 0.1-3 keV in the LECS, and #36-220 for 1.6-10 keV in the 
MECS) have been then rebinned using the grouping templates 
available at Be/?/?oSAX-SDC 10 The background has been eval- 
uated from the blank fields provided by the BeppoS AX-SDC 1 1 , 
using an extraction region similar in size and position to the 
source extraction region. The Nov. 1998 release of public cali- 
bration files, matrices and effective areas was used. 

For the PDS spectra we applied the improved screening, im- 
plementing the temperature and energy dependence of the pulse 
rise time (so-called PSA correction method). 

Due to a slight mis-match in the cross-calibration among 
the different detectors, it has been necessary to include in 
the fitting models multiplicative factors of the normalization 
(LECS/MECS and PDS/MECS ratios). The correct absolute 
flux normalization is provided by the MECS (the agreement 
between MECS units is within the 2-3% limit of the system- 
atics). The expected value of these constant factors is now 
well known and does not constitute a major/additional source 
of uncertainty: according to Fiore, Guainazzi & Grandi (1999) 
the acceptable range for LECS/MECS is 0.7-1.0, with some 
dependence on the source position in the detectors, while the 
PDS/MECS ratio can be constrained between 0.77 and 0.95. 
The latter parameter is indeed crucial in sources like Mkn 421 
where there is a significant possibility that a different compo- 
nent arises in the PDS range, whose recognition critically de- 
pends on the capability/sensitivity to reject the hypothesis that 
the PDS counts can be accounted for by the extrapolation of the 
LECS-MECS spectrum (see Section 4.3.3). 

4. 1 . Single and Broken Power Laws 

We do not discuss in detail any power law spectral models. 
For each LECS+MECS spectrum we fit the data with the single 

'A word of caution is necessary about localized features probably related 
with calibration problems yet to be solved (e.g. E~ 0.29 keV, ~ 2 keV, in 
correspondence to the Carbon edge and the Gold features due to the optics, 
respectively). 

10 The energy resolution of the LECS and MECS detectors is lower than 
the energy spacing between instrumental calibrated channels (~ 10 eV for 
LECS and 45 eV for MECS). Therefore it is compelling to rebin the spectra 
even in cases (like ours) where the statistics in each low energy channel is 
good, in order to avoid to overweight (as Xa + Xb > Xa+t,) information that 
is not truly independent. BeppoSAX-SDC templates are downloadable at 
ftp: / / www .sclc. asi.it/ pub/ sax/cal / responses /grouping/ . 
Due to the very steep spectral shape, in order 
to maintain a good statistics in each new bin, 
we altered the templates for MECS above — 7 keV, 

increasing the grouping of the original PI channels. 

"Blank fields event files were accumulated on five different 
pointings of empty fields and are available at the anonymous ftp: 
ftp : / / www .sdc. asi.it/ pub/ sax/cal/bgd/. 
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Table 1 

LECS and MECS Curved Model Spectral Fit: 1997 and 1998 
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d Since there are always three varying parameters -either two a and E% or one a, one E and E%- all the quoted errors are for 1 a for 3 parameters 
i.e. Ax 2 = 3.53), a quite conservative choice. 

b Unabsorbed Fluxes in units of 10~'° erg cnT 2 s~', in the reported energy band, whose boundaries are expressed in keV. 

c The degrees of freedom are 57 for all cases 
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Table 2 

Summary of Power Law models Results 



Model 




1997 








1998 






N H a 


x l 


xi 


cl.o.f. 


N H 11 


x 1 


X }' 


d.o.f. 


single RL. 


Galactic 


11903.1 


12.609 


944 


Galactic 


11245.0 


15.883 


708 


single RL. 


3.78±0.18 


3012.7 


3.246 


928 


3.22±0.14 


3413.2 


4.905 


696 


broken RL. 


Galactic 


1162.6 


1.275 


912 


Galactic 


997.6 


1.458 


684 


broken RL. 


2.23 ±0.46 


1083.2 


1.209 


896 


2.05 ±0.19 


866.6 


1.289 


672 


curved 


Galactic 


935.7 


1.026 


912 


Galactic 


721.1 


1.054 


684 


curved 


1.73±0.12 


917.8 


1.024 


896 


1.70 ±0.12 


705.2 


1.049 


672 



absorbing column units are 10 cm 




Energy [keV] 



FIG. 2. — Data/Model ratios for different spectral models applied to the 
1998 data. Top to bottom: single power law with (a) Galactic and (b) free 
absorbing column, broken power law with (c) Galactic and (d) free absorbing 
column, curved model, with (e) Galactic and (f) free absorbing column. The 
corresponding values of \ 2 are reported in Table 2. 



and the broken power law models 12 , both with free and fixed 13 
absorbing column density. 

Here we just show that they are an inadequate description of 
the downward curved spectra 14 . To illustrate the discrepancy, 
we took the data/model ratio for the best fit continuum model, 
for each of the individual LECS+MECS spectra (16 for 1997, 
12 for 1998), and we summed all of them, propagating the er- 

12 In Appendix C we enclose a Table reporting the best fit values for the the 
relevant parameters for the broken power law model. 

13 In the direction of Mkn 421 the Galactic equivalent absorbing column is 
N H = (1.61 ±0.1) X 10 20 cirT 2 (Lockman & Savage 1995), with the estimated 
uncertainty for high Galactic latitudes. 

14 In Fossati et al. (1998b) we presented a preliminary analysis of 1997 data, 
where LECS and MECS spectra have been fit with a variety of combinations of 
single and broken power law models to better describe the spectral curvature. 
We refer the interested reader to that work. 



rors accordingly. The results are shown in Figure 2 for the 1998 
data only, while in Table 2 we report the total y 2 value for each 
trial model, together with the (simple) average value of Nh for 
fits with free absorption. 

Neither the single nor the broken power laws provide an ad- 
equate representation of the data, not even when the value of 
the absorbing column density is left free to vary 15 -a common 
workaround used when the spectrum shows either a soft excess 
or a marked steepening toward higher energies (e.g. Takahashi 
et al. 1996). 

Indeed the description of the curvature by means of soft X- 
ray absorption is not only un-physical, but it hinders the possi- 
bility of extracting all the information from the data, by leaving 
as meaningful parameter only the higher energy spectral index 
and accounting for all the other effects by Nh- Furthermore 
this yields spectral indices that are only rough estimates, as it is 
clear that the best fit single power law does not describe the ob- 
served data in any (even narrow) energy range (see Figure 2a,b). 
Finally, in the case of Mkn 421 there is no reason to postulate 
any intrinsic absorbing component responsible for the observed 
spectral curvature (see also §4.3.2). 

4.2. The Curved Model 

Motivated by the failure of the simplest power law models, 
we developed a spectral model which is intrinsically curved 16 , 
with the aim of extracting all the information contained in the 
data, and in particular estimate the position of the peak of the 
synchrotron component, one of the crucial quantities in blazars 
modeling. 

A spectral model providing an improved description of the 
continuum shape also allows a more sensitive study of dis- 
crete spectral features such as absorption edges, whose pres- 
ence has been long searched in the soft X-ray blazars spectra 
(e.g. Canizares & Kruper 1984; Sambruna et al. 1997; Sam- 
bruna & Mushotzky 1998). 

We started from the following general description of a con- 
tinuously curved shape (see e.g. Inoue & Takahara 1996, 
Tavecchio, Maraschi & Ghisellini 1998): 



F(E) = KE~ 




(1) 



''Although the plot shown in Figure 2d could appear as a good fit, the sig- 
nificance of the small deviations above a few keV is indeed high, yielding an 
unsatisfactory fit. 

16 In Appendix A we discuss an important caveat concerning the definition of 
curved models. 
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where a_oo and a +oa are the asymptotic values of spectral in- 
dices for E <C Eg and E Eb respectively, while Eg and / 
determine the scale length of the curvature. 

We re-expressed this function in terms of the spectral indices 
at finite values of E, which characterize the local shape of the 
spectrum, instead of the asymptotic ones, which do not have 
any direct reference to the observed portion of the spectrum. 

The spectral model is then expressed in a form such that the 
available parameters are {E\ , a \ , Ei , «2 , £b , /) instead of (a-00 , 
a +00 , Eb, f) (for more details see Appendix B). As we have 
two extra parameters, for a meaningful use of this spectral de- 
scription we have to fix one for each of the pairs {E\,a\) and 
{E2,ct2). Eventually this degeneracy turns out to be a power- 
ful property of this model, because it allows us to derive the 
spectral index at selected energies (setting E\ at the preferred 
values), or even more interestingly to estimate the energy at 
which a certain spectral index is obtained (setting ot\ at the de- 
sired value). The most relevant example of this latter possibility 
is the determination of the position of the peak (as seen in vF v 
representation) of the synchrotron component Speak -if it falls 
within the observed energy band- and the estimate of the asso- 
ciated error. This can be obtained by setting one spectral index, 
i.e. ct\ = 1, and leaving the corresponding energy E\ free to vary 
in the fit: the best fit value of E(a =1) gives £ pea k- 

4.3. Results 

Most of the spectral analysis with the curved model has been 
performed by keeping the absorbing column fixed to the Galac- 
tic value, to remove the degeneracy between the effects of a 
variable absorption and intrinsic spectral curvature. We will 
briefly comment on the fits with a variable absorbing column in 
the next section. 

In order to determine the values for the really interesting pa- 
rameters, we tried a set of values for the parameter /, ranging 
from 0.5 to 3 (with a step of 0.5). All the other parameters were 
left free to vary during these step-fits. We then selected the 



value of / yielding the minimum total \ 2 - The best fit value 
for 1997 is / = 1, while for 1998 / = 2, as expected from the 
stronger curvature of the spectra. In Table 1 we report the spec- 
tral parameters, together with the 0.2-1, 2-10 and 0.1-10 keV 
fluxes (de-reddened) and \ 2 values, for each of the time sliced 
spectra. The global Data/Model ratio is shown in Figure 2e, 
and the total x 2 is reported in Table 2. The curved model is 
strongly preferred also from the purely statistical standpoint: in 
fact the model with Galactic N H has a significantly better \ 2 
than the broken power law with free Nh, although it has one 
less adjustable parameter (i.e. one more degree of freedom). 

Each spectrum has been fitted a few times in order to de- 
rive spectral indices at several energies (0.5, 1, 5, 10 keV), and 
also Ep ea k. For consistency we checked each time that not only 
the value of the \ 2 remained the same, but also all the other 
untouched parameters took the same value and confidence in- 
tervals. The fits are indeed very robust in this sense. 

The main new result is that we were able to determine the 
energy of the peak of the synchrotron component, and to assign 
an error to it. This has been possible with reasonable accuracy 
for both 1997 and 1998 spectra, with a couple of cases yielding 
only an upper limit for the peak energy 17 . 

An example of the remarkable spectral evolution during the 
1998 flare is shown by the deconvolved spectra in Figure 3, 
which also illustrates the strong convexity of the spectrum and 
the well defined peak. 

In 1997 the source was in a lower brightness state, with an 

average X-ray spectrum softer (at all energies Aagj — 0.4), 

and a peak energy 0.5 keV lower. 

More globally, the analysis shows that there is a clear rela- 
tion between the flux variability and the spectral parameters, 
both in 1998 and, albeit less strikingly, in 1997. Particularly 
important is the correlation between changes in the brightness 
(even the small ones) and shifts of the peak position, as the latter 
carries direct information on the source physical properties. In 
Figure 4a,b, a at 5 keV and E pea k for 1997 and 1998 are plotted 
versus the 0.1-10keV flux. The source reveals a coherent spec- 
tral behavior between 1997 and 1998 and through a large flux 
variability (a factor 5 in the 0.1-10 keV band). Both the peak 
energies and the spectral indices show a tight relation with flux, 
the latter being described by E pea k oc F e , with e = 0.55 ± 0.05. 

4.3.1. Hard Lag in 1998 spectra 

One further important finding of the time resolved spectral 
study is the signature of the hard lag which strengthens what 
already inferred from the timing analysis (see Paper I). In fact 
the local spectral slope at different energies indicates that the 
flare starts in soft X-rays and then extends to higher energies. 

A blow up of the 1998 flare interval is shown in Figure 5: 
spectral indices at 1 and 5 keV, and E pea k are plotted versus time 
in the left hand plots, and versus flux in the right hand ones. The 
synchrotron peak shifts toward higher energies during the rise 
and then decreases as soon as the flare is over. The spectral 
index at 1 keV reflects exactly the same behavior, as expected 
being computed at the energy around which the peak is moving 
(and in fact it moves in a narrow range around a = 1). On the 
contrary, the spectral shape at 5 ke V does not vary until a few ks 
after the peak, and only then -while the flux is decaying and the 
peak is already receding- there is a response with a significant 

"These are the cases "35-37" and "41^12" of 1997. Actually the formal 
fitting yields a confidence interval, but with the lower extreme falling outside 
of the observed energy range. 



6 



Evolution of Mkn 421: II. Spectral Analysis and Phsyical Constraints 



e 2 




3 4 5 6789 10 

Fo.i-iouev [1°- 10 erg/cmVs] 



> 

CD 




3 4 5 6789 10 

F .i-ioi»v [10- 10 erg/cmVs] 



FIG. 4. — (a) spectral index at 5 keV and (b) synchrotron peak energy plotted 
versus the de-absorbed 0.1-10 keV flux, both for 1997 and 1998 datasets. The 
dashed line across the diagram (b) represents the best fitting power law, with a 
slope e = 0.55. 



hardening of the spectrum. The right side diagrams illustrate 
this behavior directly in terms of flux: the evolution of the peak 
energy and the 1 keV spectral shape follow closely that of the 
flare, while at 5 keV the spectrum does not change shape until 
after the peak. 

Thus the spectral evolution at higher energies develops dur- 
ing the decay phase of the flare, tracing a loop in the a vs. Flux 
diagram, with a significant hardening of the spectrum. Due to 
the presence of the hard lag, the loop is traced in the opposite 
way with respect the evolution that is commonly observed dur- 
ing flares in HBL, that is a spectral hardening during the rise 
of the outburst, followed by a softening during the decay (see 
however Sembay et al. 1993, Catanese & Sambruna 2000). 

This behavior -usually accompanied by a soft lag- produces 
a clockwise loop in a vs. Flux. 

4.3.2. Absorption Features: Nu and edges 

We also fit the 1997 and 1998 spectra with the curved spectral 
model leaving the N H free to vary. The only constraint that 
we imposed is on the exponent /, which is set at the best fit 
values of / = 1 and / = 2 for 1997 and 1998, respectively. The 
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improvement of the global fit 18 is not statistically significant. 
The average best fit absorbing column is indeed very close to 
and consistent with the Galactic one, leading us to conclude that 
there is no requirement whatsoever for an additional absorbing 
column to model the data. 

We searched for the presence of absorption edges at low ener- 
gies (E £ 1 keV), both in 1997 and 1998 data. We added to the 
curved continuum an absorption edge, with a first guess fit at 
E~ 0.4-0.5 keV, as even in the case of the best curved model 
with free N H there are residues around this energy (see Fig- 
ure 2f). The most delicate issue is the calibration problem re- 
lated with the instrumental Carbon absorption edge at 0.29 keV. 
Even the smallest uncertainties (~ 1-2%) in the knowledge of 
the instrument response around this energy can induce signif- 
icant spurious features, typically at slightly higher energy (see 
Fossati & Haardt 1997). 

There is not convincing evidence for discrete absorption fea- 
tures. In 5/28 cases the fit routine has not found an edge, while 
there are indeed a few spectra (8/28) in which the fit definitely 
put the edge at 0.27-0.30 keV, to account for systematic devia- 
tion due to calibration uncertainties. There are then 15/28 cases 
where the edge energy is around 0.4-0.6 keV, but again only 
in a handful of cases the best fit energy E e d ge is significantly 
different from 0.29 keV. 

Summarizing, there are only three individual spectra (out 
of 28) for which the inclusion of the edge is formally statisti- 
cally significant (> 95%), according to the F-test. These three 
cases do not stand out in the sample for any other property, 
such as brightness or curvature of the continuum. Best fit en- 
ergies and optical depths are typically E e d ge — 500+fg eV, and 
t ~ 0.4^o 2 _0 ' 6 . However, considering that the sample comprises 
28 spectra, with the significance threshold set at 95%, three pos- 

18 Ax 2 = 17.9 for 16 additional d.o.f. for 1997 and Ax 2 = 15.9 for 12 addi- 
tional d.o.f. for 1998, corresponding to an average A\ 2 — 1.2 per spectrum for 
a change in d.o.f. from 57 to 56. 
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itive detections can not be regarded as a compelling evidence 
for the presence of absorption edges. 

4.3.3. PDSdata 

Covering the range above ~ 12 keV, the PDS instrument 
could detect the IC component, which might start to dominate 
the emission in this band. It has always proven to be very diffi- 
cult to detect/constrain this component, although it is expected 
to have a very hard spectrum, that should make it easier to dis- 
entangle from the very steep tail of the synchrotron component. 
There are a few objects, the so-called intermediate BL Lacs, 
in which the cross-over between synchrotron and IC occurs in 
the 0.1-10 keV band, for which it has then been possible to 
clearly detect the hard IC power law, with typical a ~ 0.3-0.7 
(e.g. Tagliaferri et al. 2000 for ON 231; Giommi et al. 1999 
for S5 0716+714; Sambruna et al. 1999 for BL Lac). The 
good sensitivity of the PDS yielded a convincing evidence of IC 
emission in the case of PKS 2155-304 (Giommi et al. 1998). 
The preliminary analysis of the Mkn 421 1997 data (Fossati et 
al. 1998b) suggested the presence of a hard component. We 
therefore focused on the search for the IC component for both 
of 1997 and 1998 PDS datasets. 

For 1998, when the source was brighter, we accumulated 
PDS spectra according to the same partition of the light curves 
used for the 0.1-10 keV data. In most of the individual spec- 
tra there is a significant (>3-cr) detection only up to 40 keV. 
The light curve for the 12-26 keV band is shown in Figure 3 of 
Paper I. 

We then restricted the analysis to only one spectrum for 1997, 
integrating over the whole campaign, and two for 1998, one for 
April 21 st , and one for April 23 rd . The integration times are 
54.8, 21.1 and 25.2 ks respectively (see Paper I). We grouped 
the data in 6 channels spectra with boundaries at 12, 18, 27, 40, 
60, 90 and 130 keV. 

In the 1997 spectrum there is a positive detection up to 
90 keV. For April 21 st -corresponding to the synchrotron flare- 
there is a strong signal up to 60 keV, while, surprisingly, in the 
April 23 rd spectrum there is a ^3-cr detection in each individual 
bin up to 130 keV. From 12 to 40 keV the April 21 st count rates 
are significantly higher than those of April 23 rd , but while the 
former continue to steeply decrease with increasing energy, the 
latter flatten above 40 keV, with a cross over in the 40-60 keV 
bin. As the integration times for these two spectra are similar 
(and in fact statistical uncertainties on the count rate are of the 
same order), this cannot account for this difference 19 . 

The results of single power law fits are reported in Table 3: 
clearly the April 23 rd spectrum is significantly harder. More- 
over, the spectral indices for the first two cases are consistent 
with those derived from the LECS and MECS datasets and the 
general picture of a continuously steepening spectrum, while 
this does not hold for the third case. 

While this is already interesting, there is one additional, 
somehow unexpected, finding: apparently the April 23 rd harder 
spectrum cannot be simply interpreted as the detection of the 
IC component when the synchrotron tail recedes. In fact a hard 
power law at this flux level -if present- would have been de- 

19 Also, we are not aware of any intra-observation "background fluctuations" 
to be taken into account in the standard PDS data reduction (which has been 
performed by following carefully the prescriptions given by the instrument 
team, for all the details refer to Chiappetti et al. 1999), nor of any long term 
"background fluctuations" that could be responsible for the observed variation. 
No such thing was reported neither by the PDS instrument team, nor by the 
BeppoSAX Science Data Center. 



tected also on April 21 st , while the count rates in the higher 
energy bins for April 21 st set a tight upper limit. This behav- 
ior might be ascribed to either a different independent origin 
or a delayed response of the IC component with respect to the 
synchrotron one (indeed some evidence of a possible delayed 
decay of the PDS with respect to the MECS emission can be 
seen in the light curves). 

Finally we tried to quantify the presence of the IC power law 
in the three datasets by fitting MECS and PDS data together. 
We proceeded in the following way: 

• We accumulated LECS and MECS spectra over the same 
intervals of the PDS ones, i.e. one cumulative spectrum 
for each day of 1998, and one for the whole 1997. 

• We fit the LECS+MECS spectra in order to constrain the 
average synchrotron spectrum. 

• We then fixed all the parameters of the curved model to 
the LECS+MECS best fit values, and added to the model 
a hard power law component 20 with fixed spectral in- 
dex (we tried the values aic = 0.4,0.5,0.6,0.7). We also 
added an exponential cut-off to the synchrotron compo- 
nent setting its two parameters, i.e. cut-off and e-folding 
energies, to 10 and 40 times the best fit value of E pea k. 

• We used this partially constrained model to fit the MECS 
(above E> 5 keV) and PDS spectra. The two 1998 
datasets have been fit jointly with PDS/MECS intercali- 
bration factor free to vary, but set to be the same for both, 
and of course independent normalization. 

The results are reported in Table 4, for the arc = 0.5 case (the 
values relative to the other aic are not significantly different). 
The best fit PDS/MECS normalization ratio is 0.88. 

There seems to be a significant change in the flux level of the 
expected IC component between 1998 April 21 st and 23 rd . The 
best fit value for the normalization for the first day is zero, while 
for the second day a non-zero IC flux is definitely required. In 
order to test the likelihood of the apparent variation, we fit the 
data constraining the normalization to be the same for the two 
datasets (third case in Table 4) and the F-test probability to 
obtain by chance the observed change in \ 2 is ~ 0.01 (A^ 2 = 
6.08, d.o.f. from 30 to 31). 

The basic result for 1997 is that the IC component is defi- 
nitely detected. Its brightness level is comparable with the av- 
erage of 1998, while the synchrotron component shows a sig- 
nificant variation. 

It should be noted that the synchrotron fluxes reported in Ta- 
ble 4 are not the best reference because at those energies, well 
above the synchrotron peak, the spectrum is very steep and even 
the smallest shift of the peak energy translates in a big change in 
the flux. Nevertheless even when considering the energy range 
0. 1-10 keV, which includes the synchrotron peak, the variation 
between 1997 and 1998 is of the order of a factor of 2-3 (e.g. 
Figure 5). The sampling of the 1997 light curve is not good 
enough to enable us to investigate any possible variation of the 
PDS spectrum with the source brightness. 

20 We used the pegpwrlw (power law with pegged normalization) model 
of XSPEC that provides a robust measure of the normalization. In fact this is 
defined over a selectable energy band, instead of the monochromatic value at 
1 keV that can be strongly correlated with other parameters if 1 keV does not 
fall close to the logarithmic median of the band spanned by the data. 
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Table 3 

pds single power law fits: 1997 and 1998 



Count rates 



Dataset 


Band 


Exposure 
(ks) 


12^10 keV 

(cts/s) 


40-90 keV 
(cts/s) 


a" 


FlkeV 

0*Jy) 


x 2 


cLo.f. 


1997 


12-90 keV 


54.8 


0.177 ± 0.021 


0.045 ± 0.019 


9 17+0.48 
Z - -0.41 


1.52 


2.96 


3 


1998 April 21 st 


12-90 keV 


21.1 


0.698 ± 0.028 


0.039 ± 0.023 


9 9-1+0.14 
Z - ZJ -0.13 


5.76 


1.76 


3 


1998 April 23 rd 


12-90 keV 


25.2 


0.293 ± 0.025 


0.081 ±0.021 


i tc+0.26 
-0.24 


1.35 


6.95 


4 




12-130 keV 








| 54+0.28 
1J -0.25 


1.54 


2.95 


3 



a Errors are for 1 a for 1 parameter (i.e. A% 2 = 1.0). 



Table 4 

High energy component constraints from MECS+PDS spectral fits 





High Energy Component Flux a ' b 


Synchrotron Flux a ' c 






Dataset 


12-40 keV 


40-100 keV 


12-40 keV 


40-100 keV 


x 2 


cLo.f. 


1997 




12.8+i;? 


7.8 


0.2 


17.92 


15 


1998 April 21 s ' 


< 13.0 d 


< 14.9 d 


78.4 


8.6 


28.71 


30 


1998 April 23 rd 


13.6+33 


15.8+11 


24.0 


1.2 






1998 April 21 s ' + 23 rd 


10.1+3 1 


11.8+" 


78.4/24.0 


8.6/1.2 


34.79 


31 



d Fluxes are in units of 10~ 12 erg cirT 2 s -1 . Errors are for 1 a for 1 parameter (i.e. A\ 2 = 1.0). 

''Flux computed by the best fit model in the hard power law. 

c Flux computed by the best fit model in the curved model component. 

d 90% upper limit. 
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5. DISCUSSION 

Before focusing on the modeling of the temporal and spectral 
behavior illustrated in this paper and in Paper I, let us consider 
a few issues which we consider of particular relevance for the 
interpretation of the origin of blazar variability. 

5.1. Spectral Variability, Synchrotron Peak and IC component 

Soft and hard X-ray bands show a different behavior. This 
might be attributed to the contribution from both the syn- 
chrotron and inverse Compton components at these energies. 
As electrons emitting X-rays through the two processes would 
have different energies, the two components are not expected to 
vary on the same timescales and in phase. Indeed on one hand 
the constraints provided by the analysis of the PDS data suggest 
that the synchrotron component constitutes the dominant con- 
tribution to the flux in the 12-26 keV band on 1998 April 21 st 
(during the flare) with a light curve somewhat different from 
that of the softer synchrotron X-rays. On the other hand, the 
evidence supporting the possibility of a substantial variation in 
the IC component during the 1998 observations is extremely in- 
teresting. We can rule out the possibility that the flux level of 
the hard component is approximately constant between the two 
observations of 1998 (with a 99% confidence), and thus that its 
detectability only depends on the the brightness level of the tail 
of the synchrotron emission. 

5.2. The Flaring + Steady Components Hypothesis 

A further interesting issue raised in Paper I in connection 
with the temporal analysis is the presence and role of quasi- 
stationary emission. Our results support the view that the 
short-timescale, large-amplitude variability events could be at- 
tributed to the development of new individual flaring compo- 
nents, giving rise to a spectrum out-shining a more slowly vari- 
able contribution. 

Indeed there is growing evidence that the overall blazar emis- 
sion comprises a component that is (possibly) changing only on 
long timescales and that does not take part in the flare events 
(see for instance Mkn 501, Pian et al. 1998, and Pian et al. 
1999 for the case of the UV variability of 3C 279). 

The deconvolution of the SED into different contributions 
would allow the study of the evolution of the flaring compo- 
nent possibly clarifying the nature of the dissipation events oc- 
curring in relativistic jets -in particular the particle accelera- 
tion mechanism- and the modality and temporal characteristics 
of the initial release of plasma and energy through these colli- 
mated structures. 

The deconvolution is however still difficult to achieve with 
the available temporal and spectral information. In paper I we 
discussed the results of our attempt to measure directly the rela- 
tive contributions of what we identify as steady and flaring com- 
ponents (parameter 1Z), from the characteristics of the flare de- 
cay. The result is very interesting, although it does not provide 
a useful handle on the evolving spectral properties of the flare 
itself. A more feasible approach is to assume specific forms of 
the flare evolution and test them against the data. One of the 
simplest working hypothesis is to assume that the flare evolu- 
tion can be reproduced by the time dependent flux and energy 
shifts of a peaked spectral component with fixed spectral shape 
(e.g. Krawczynski et al. 2000). In this context it is relevant 
to remind that the maximum of the synchrotron emission cor- 
responds to the energy at which most of the energy in particles 
is located. Its value and temporal behavior thus trace the evo- 



lution of the bulk of the energy deposition into particles by the 
dissipation mechanism. 

It is therefore possible that the variability characteristics 
(both temporal and spectral) might chiefly depend on the (par- 
ticles/photons) energy relative to the synchrotron peak. In 
particular, the results on T s hort (Paper I) can be either related 
to the dominance of the light crossing time over the cooling 
timescales, or alternatively to the different position of the sam- 
pled energies with respect to the synchrotron peak as above 
the peak the timescales are shorter. In any case the fact that 
timescales are similar at the different energies possibly indi- 
cates the evolution of a (observed) fixed-shape flaring compo- 
nent. 

5.3. Sign of the lag 

One of the main new results of the temporal analysis pre- 
sented in Paper I is the significant detection of a hard lag, oppo- 
site to the behavior commonly detected in HBL. Interestingly, 
a recent study by Zhang et al. (1999) on the source PKS 2155- 
304 has shown the presence of an apparent inverse trend be- 
tween the source brightness level and the duration of the (soft) 
lag. And intriguingly, the 1998 flare represents the bright- 
est state ever observed from Mkn 421, thus possibly suggest- 
ing that an extrapolation of the above phenomenological trend 
might even account for a hard lag. This possibility is currently 
under study (Zhang et al., in preparation). 

Again the relative contribution of two components might be 
responsible for this trend. During the most intense events the 
flaring component would completely dominate over the quasi- 
stationary one, progressively shifting to higher energies (hard 
lag), while in weaker flares the varying component would ex- 
ceed the steady one only at energies higher than the peak, where 
the spectrum steepens, producing an observed soft lag behavior. 

However, one should also consider that in different bright- 
ness states the different observed energies correspond to differ- 
ent positions with respect to E pea k (e-g- comparing the 1997 and 
1998 datasets we showed that its value was typically a factor 2 
different), requiring a more subtle analysis before deriving in- 
ferences on the lag-brightness relationship. 

5.4. Physical Interpretation: the signature of particle 
acceleration 

Let us now focus on the interpretation of the two main and 
robust results of this work, namely the hard lag and the evolu- 
tion of the synchrotron peak. 

The occurrence of the peak at different times for different en- 
ergies is most likely related to the particle acceleration/heating 
process. Although models to reproduce the temporal evolution 
of a spectral distribution have been developed (e.g. Chiaberge 
& Ghisellini 1999, Georganopoulos & Marscher 1996; Kirk, 
Rieger & Mastichiadis 1998), they mostly do not consider the 
role of particle acceleration (see however Kirk et al. 1998). 

In order to account for the above results, we thus introduced 
a (parametric) acceleration term in the particle kinetic equation 
within of the time dependent model studied by Chiaberge & 
Ghisellini (1999). In particular, this model takes into account 
the cooling and escape in the evolution of the particle distribu- 
tion and the role of delays in the received photons due to the 
light crossing time of different parts of the emitting region. 

The model includes the presence of a quiescent spectrum, 
which is assumed to be represented by and thus fitted to the 
broad band spectral distribution observed in 1994 (Macomb et 
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al. 1995). To this a flaring component is added and this is con- 
strained by the observed spectral and temporal evolution (the 
parameters for both the stationary and variable spectra are re- 
ported in Table 5). The values of the model parameters are 
very close to those derived in Maraschi et al. (1999) in a some- 
what independent way, i.e. trying only to fulfill the constraints 
provided by the X-ray and Te V spectra, simultaneous but aver- 
aged over the whole flare. On this basis it is very likely that the 
model discussed here satisfies the TeV constraints, which are 
however not directly taken into account (the computation of the 
TeV light curve in this complex model would require to treat 
the non-locality of the IC process, which is beyond the scope 
of this paper). 

Clearly a parametric prescription does not reproduce a priori 
a specific acceleration process, but we rather tried to constrain 
its form from the observed evolution. The main constraints on 
the form of the acceleration term are the following: 

a) Particles have to be injected at progressively higher ener- 
gies on the flare rise timescale to produce the hard lag. 

b) Globally the range of energies over which the injection 
occurs has to be narrow, to give rise to a peaked spectral 
component. 

c) A quasi-monochromatic injection function for the flux at 
the highest energies to reach its maximum after that at 
lower ones (as opposed to e.g. a power law with increas- 
ing maximum electron Lorentz factor). 

d) The emission in the LECS band from the particles which 
have been accelerated to the highest energies (i.e. those 
radiating initially in the MECS band) should not exceed 
that from the lower energy ones, as after the peak no fur- 
ther increase of the (LECS) flux is observed; this in par- 
ticular requires for the injection to stop after reaching the 
highest energies. 

e) The total decay timescale might be dominated by the 
achromatic crossing time effects, although the initial fad- 
ing might be determined by different cooling timescales. 

The acceleration term has then been described as a Gaussian 
distribution in energy, centered at a typical particle Lorentz fac- 
tor 7 c (t) and with width a = 0.0l7 c (t), exponentially evolving 
with time: 7 c (t) oc e~ {lm "~ l \ where t max corresponds to the end 
of the acceleration phase, when the maximum 7 c , max is reached. 
The injected luminosity is assumed to be constant in time. 

The duration of the injection, assumed to correspond to the 
light crossing time of the emitting region, is such that it mim- 
ics the passage of a shock front (i.e. a moving surface passing 
through the region in the same time interval, see Chiaberge & 
Ghisellini 1999). 

It should be also noted that, within this scenario, the sym- 
metry between the rise and decay in the soft energy light curve 
seems to suggest that if the energy where most of the power is 
released is determined by the balance between the acceleration 
and cooling rates, at this very same energy the latter timescales 
are comparable to the light crossing time of the region. 

The predictions of the model are shown in Figures 6 and 7, in 
the form of light curves and spectra at different times, respec- 
tively. In particular, in Figure 6 the light curves (normalized 
over the stationary component) at the two centroid energies of 
the soft and hard bands considered in the data analysis are re- 
ported. The presence of a hard lag is clearly visible as well as 
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FIG. 6. — Light curves computed from the model, at the centroid energies 
of the soft and hard band considered in the data analysis. The vertical axis 
represents the variations normalized to the flux of the stationary component. 
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the larger variability amplitude at higher frequencies. Also the 
spectral evolution, reported in Figure 7 in vF v in the observed 
energy range, seems to be at least in qualitative agreement with 
what observed (compare with Fig. 3). 

A posteriori it is important to notice that the inclusion of a 
quasi-stationary component and of high energy emission (from 
synchrotron self-Compton) are indeed crucial in order to repro- 
duce the retarded spectral variability at energies above a few 
keV. This constitute a further difference with respect to the dis- 
cussion of Maraschi et al. (1999), where a simplified descrip- 
tion of the time decay was used, based on evidence for an en- 
ergy dependence of the flare decay timescale (obtained model- 
ing the decay with no baseline). 

6. CONCLUSIONS 



Fossati et al. 



11 



Table 5 
Model Parameters 



Parameter 


Parameter Description 


stationary 


flaring 


R 


source dimension (cm) 


2 -10 16 


2-10 16 


B 


magnetic field (G) 


0.12 


0.03 


tesc 


escape timescale (R/c) 


6 


1.5 


f- - 
>nj 


compactness injected as particle energy 


6-irr 5 


3 -10- 4 


S 


slope of the injected stationary power law 


1.6 




7max 


maximum Lorentz factor of the injected stationary power law 


3 10 5 




O'cmax 


maximum Lorentz factor of Gaussian injection for the flaring component 




7 -10 s 


8 


beaming factor 


18 


18 



BeppoSAX has observed Mkn 421 in 1997 and 1998. We 
analyzed and interpreted the combined spectral and temporal 
evolution in the X-ray range. During these observations the 
source has shown a large variety of behaviors, both concerning 
the X-ray band itself and its variability properties with respect 
to the 7-ray one, providing us with a great wealth of infor- 
mation, but at the same time revealing a richer than expected 
phenomenology. 

Several important results follow from this work: 

(a) The X-ray and 2 TeV light curves peak simultaneously 
within one hour, although the halving time in the TeV 
band seems shorter than those at LECS and MECS ener- 
gies (see Maraschi et al. 1999). 

(b) The detailed comparison of 0.1-1.5 keV and 3.5-10 keV 
band light curves shows that the higher energy band lags 
the softer one, with a delay of the order of 2-3 ks. This 
finding is opposite to what has been commonly observed 
in HBL X-ray spectra (see Paper I). 

(c) Moreover, extracting LECS+MECS spectra for <^ 5 ks in- 
tervals, we were able to follow in detail the spectral evo- 
lution during the flare. For the first time we could quan- 
titatively track the shift of the peak of the synchrotron 
component moving to higher energy during the rising 
phase of the flare, and then receding. 

(d) An energy dependence of the shape of the light curve dur- 
ing the flare has been revealed: at low energies the shape 
is consistent with being symmetric, while at higher ener- 
gies is clearly asymmetric (faster rise) (see Paper I). 

(e) Evidence has been found for the presence of the IC com- 
ponent, and more importantly for its substantial variabil- 
ity, which is possibly delayed with respect to the syn- 
chrotron one. 

These findings provide several temporal and spectral con- 
straints on any model. In particular, they seem to reveal the 
first direct signature of the ongoing particle acceleration mech- 
anism, progressively "pumping" electrons from lower to higher 
energies. The measure of the delay between the peaks of the 
light curves at the corresponding emitted frequencies thus pro- 
vides a tight constraint on the timescale of the acceleration pro- 
cess. 

Indeed, within a single emission region scenario, we have 
been able to reproduce the sign and amount of lag by postulat- 
ing that particle acceleration follows a simple exponential law 



in time, stops at the highest particle energies, and lasts for an 
interval comparable to the light crossing time of the emitting 
region. If this timescale is intrinsically linked to the typical 
source size, we indeed expect the observed light curve to be 
symmetric at the energies where the bulk of power is concen- 
trated and an almost achromatic decay. The same model can ac- 
count for the spectral evolution (shift of the synchrotron peak) 
during the flare. 

The other very important clue derived from the analysis is 
the presence of a quasi-stationary contribution to the emission, 
which seems to be dominated by a highly variable peaked spec- 
trum, possibly maintaining a quasi-rigid shape during flares. 
The decomposition of the observed spectrum into these two 
components might allow us to determine the nature and modal- 
ity of the energy dissipation in relativistic jets. 

We are grateful to the BeppoSAX Science Data Center (SDC) 
for their invaluable work and for providing standardized prod- 
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edge the Italian MURST for financial support. This research 
was supported in part by the National Science Foundation un- 
der Grant No. PHY94-07194 (AC). Finally, GF thanks Cecilia 
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APPENDIX 

A. A CAVEAT ON CONTINUOUSLY CURVED SPECTRAL 
MODELS 

The increasing quality of the available X-ray spectral data, 
both in terms of signal-to-noise ratio and energy resolution, has 
made necessary to consider more complex fitting models. One 
of the most interesting continuum feature in the X-ray spec- 
tra of blazars is the curvature of the synchrotron component, as 
good quality data could enable us to estimate the energy of the 
emission peak. As clearly the traditional single or broken power 
law models are unsatisfactory, we developed the curved spectral 
model presented in this paper: a simple analytic expression rep- 
resenting a continuous curvature, two "pivoting" points -useful 
for analysis purposes- asymptotically joining two power law 
branches, whose slope is completely determined by the behav- 
ior of the function at the pivoting energies. 

Another increasingly popular model has been introduced 
by Giommi et al. (1998) (to reproduce BeppoSAX data of 
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PKS 2155-304, and then used also for Mkn 421 by Guainazzi 
et al. 1999 and Malizia et al. 1999), but unfortunately this is 
affected by a problem that can lead to misleading results. 

This more generally occurs for all curved continuum models 
defined as: 

F(E) = E~ &(E) (Al) 

where the function a(E) is suitably defined to smoothly change 
between two asymptotic values a\ and a 2 . The specific choice 
by Giommi et al. (1998) is: 



a, = 0.75, a,= 1.5, E„=3.5 keV, fi=l 



a(E) = f(E)a l + (l-f(E))a 2 



with 



f(E) 



(A2) 



(A3) 



Thus a — > ct\ for E <C Eq, and a — > a 2 for E ^> Eq. While this 
exactly the behavior of the spectral index function, this does 
not correspond to the real spectral index of the function F(E), 
conventionally defined as: 



a (E) = dlnF(E ^ 
dlnE 

In fact for F(E) = £- a(£) this yields 
, da(E) 



a(E) = a(E) + InE 



dlnE 



&(E) + 5 & (E) 



(A4) 



(A5) 



thus including an unavoidable term 5a(E), present because of 
the not vanishing (by definition) derivative of a(E). 
For the above specific choice of a(E): 

5 & {E) = (3 —InE (ai-a 2 )e- E/Ea ( 1 -e" £ M & ' . (A6) 

Eq \ ) 

The real spectral index then takes a value (sensibly) different 
from the one expected over a large range of energies. 

The effect can be better illustrated with an example. Let us 
consider a set of spectral parameters that would well reproduce 
the spectrum of the peak of the 1998 flare of Mkn 421: ol\ = 
0.75, a 2 = 1. 5,£o = 3.5 keV, (3 = 1.0. In Figure A8 we show the 
spectrum together with a(E) and a(E) as derived by the model 
over an energy interval broader than that effectively covered by 
data. For reference we also plotted the values of a derived from 
our curved model fits, to show that the real spectral index a(E) 
matches them. A few features emerge: 

1. according to a(E) between Eq and 10 keV there is an 
apparent wiggle in the spectrum, which steepens beyond 
the supposedly asymptotic value a 2 before reaching it at 
about 30 keV. 

2. within the 0.1-10 keV range a(E) and a(E) never match 
each other except for E = 1 keV, when 5a(E) vanishes. 

3. the range spanned by the true spectral index is broader, 
and in particular it is always steeper than the value ex- 
pected on the basis of a(E) in the whole band 1-20 keV. 
Aa can be as much as +0.42 (at E ~ 5.5 keV). 

4. at no energy in the 0.1-10 keV band the actual spectrum 
has a slope equal to either ct\ or a 2 , and thus these two 
parameters do not have any real descriptive meaning (this 
is strikingly illustrated by the displacement of the real 
data points from the expected curve). 
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FIG. A8. — Example of the fitting spectral model used by Giommi et 
al. (1998), for a choice of model parameters suited to reproduce the spectral 
properties of the 1998 flare (shown on top of the figure box). The solid line is 
the spectrum (scale on the left Y-axis). The dot-dashed line is the predicted 
spectral index a(E), and the dashed line is the "true" spectral index of the re- 
sulting F(E) (both with the scale on the right Y-axis). The squares represent 
the values of the spectral index at the same energies reported in Table 1 . Verti- 
cal bars indicate a few relevant reference energies, namely 0.1 and 10 keV, that 
are the limits of the band covered by LECS and MECSs, 1 keV as the energy at 
which the model is normalized to unity, and Eq. 



5. the total change of slope in the 0. 1-10 keV interval would 
be expected Aa ~ 0.68 while it actually turns out to 
be 1.04 (and even more dramatically between 0.5 and 
5.0 keV, Aa ~ 0.46 while for the data it is 0.86; see Ta- 
ble 1). 

Summarizing, there is no way to define a curved continuum 
model by means of a suitable function describing the power 
law exponent. 

B. MORE DETAILS ON THE ADOPTED CURVED MODEL 
The relationship between a(E) and a-oo and a +00 is: 



a(E) = - 



yielding: 



dlnF(E) 
dln(E) 



f 



(E/Eb) 
l + (E/E B ) J 



a +0 



B(\+A)ai-A(\+B) a 2 
(1 + 5) a 2 -(l+A)a 1 



(Bl) 

(B2) 
(B3) 



where A = (Ei/E B ) J , andB = (E 2 /E B ) } 

C. MORE DETAILS ON THE BROKEN POWER LAW FITS. 
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Table C6 

LECS and MECS Broken Power Law Model Spectral Fit: 1997 and 1998 



Orbits Galactic absorbing column" free absorbing column b 





Ql 


Q2 


Ebreak 

(keV) 


x 2 , 


N H 
(10 20 cm" 2 ) 


"1 


Q2 


Ebreak 

(keV) 


xl 


1997 April 29 th to May 5 th 


1-3 
4-7 


96 +<u)i 
u - -0.07 

1 01+ 004 


j 7J+0.04 

i 79 +o:o3 

1 -'^-0.03 


1 77+O.I8 
1 - J '-Q.19 
1 23+ 012 


1.71 
1.55 


2 50+ - 32 
z - JU -0.30 

1 96+ 037 


1 26 +(M)8 
1 - zo -0.08 
1 14+0.13 


1 77+O.O6 
1 -"-0.05 
1 on+0.04 
LSU -0.04 


9 95+O.23 
z - z -0.23 
1 70+O.27 


1.37 
1.51 


8-10 


93+0-07 
y -0.09 


1 77+0.04 
1 -"-0.04 


1 ig+0.17 


1.38 


9 99+0.43 
z - yz -0.43 


1 40+0- 10 


1 95+0- 13 


2 90+0- 48 
z - w -0.51 


1.23 


11-12 


1 04+ 005 


i 77+0.04 
1 -"-0.04 


1 41+O.I8 


1.26 


1 Q7+0.45 


1 le 4 *' 14 

1 - 1D -0.16 


1 77+O.O5 

^"-o.os 


1 57+O.44 
-0.24 


1.22 


13-15 


97+0.04 


i O9+0.03 


1 26 +ou 
1, ' 4o -0.10 


1.24 


1 94+0.43 


1 09 +0 15 

1 ' u ^-0.17 


1 09+0.04 

1 - SZ -0.04 


1 35+ ' 24 

1,J -0.17 


1.21 


16-18 


1 00 +0 07 


^Tol 


1 14+O.n 


1.10 


1 62" H)68 

1 - OZ -0.58 


1 01+ - 25 

l.ui_ 025 


1 76+0-05 
1 -'°-0.05 


1 14+0.31 


1.12 


19-20 


Q7+0-05 
u - y/ -0.06 


1-68^ 


1 15 +(m 
1,1J -o.io 


1.41 


1 co+0.49 
O8 -0.46 


96+° 18 
u ' yo -0.20 


1 fio+0.04 
1.68_ aM 


1 14+O.I 7 
la ^-0.13 


1.44 


21-23 


94+0-04 

u -^-0.04 


i 77+O.O3 
-0.03 


1 90+0.09 

1 -^°-0.08 


1.84 


1 89 +a34 


1 04 +0 12 


1 77+O.O3 
1 - ' -0.03 


1 74+O-I6 


1.81 


24-25 


95+0-06 

u - -0.07 


1 70+0- 04 

1 -' u -0.04 


1 31 rt ' 17 


1.13 


941+0.34 

z -^-0.33 


1 95+O.O8 
i ' ZJ -0.08 


1 77+O.O7 
'■■"-0.06 


231+0.31 

z - J1 -0.28 


0.98 


26-27 


97+0.O6 


1 76 +fl04 

1 - ' -0.03 


1 9O+0.17 


1.14 


9 79+O.32 
J -0.33 


1 97+O.O7 
1 - z '-0.09 


1 on+0.08 


2 30 +0 31 
Z - J -0.28 


0.99 


28-30 


I 09+0.06 


l 77+0.04 
1 -"-0.04 


1 77+O.26 
-0.16 


1.24 


2 06+° 50 
Z ' UD -0.57 


1 25+ 018 
1,z -0.22 


1 77+O.O5 

^"-o.os 


1 54+ a92 
1 - J -0.35 


1.21 


31-32 


1 18+° 04 
lao -0.05 


2 04+° 09 
z - -0.09 


9 34+0.19 
-0.27 


0.94 


2.20^- 4 5 


1 32+ au 




9 47+O.25 
z - HJ -0.28 


0.80 


33-34 


i fiR+0 09 


2 01+° 08 

z - ul -0.08 


1 14+O.I6 
1-1 -0.13 


1.36 


3.<*£* 


1 64+°. 13 


2.28^ : ll 


7 19+O.6I 
J --0.80 


1.30 


35-37 


i 99+O.O6 
1,zz -0.07 


1 qq+0.05 


1 9^+0.17 
'•■"-0.15 


1.20 


2 56 +0 ' 38 

Z - J -0.38 


1 iro+0.10 

i - 58 -0.10 


2 15+ 013 
z>1 -0.12 


2 64 +0.39 
Z - O4 -0.50 


1.06 


38-40 


1 07+° 06 

-0.07 


1 90+0 05 


1 11+ - 11 
1-11 -0.11 


1.34 


9 09+O.34 
2.82_ Q37 


1 55+0.08 

-0.10 


2 06+°. 12 
Z - UD -0.14 




1.15 


41^12 


l 1 1+0.07 
1,1 -0.08 


1 85+° 06 

-0.06 


1 95+0.22 
' ,iJ -0.19 


0.74 


1 71+0.62 
-0.61 


1 14+O- 22 
1-1 -0.27 


1 oc+0.07 

^"-o.oe 




0.75 


1998 April 21 s ' 


1-2 
3 




1 57+ 005 

■ JJ -0.05 


1 97+O.I4 
-0.23 


0.99 
1.07 


2 01+ 66 

IT ^ 
1 - /1 -0.54 


099+0.24 

094+8+? 
u y -0.22 


1 S7+O.O6 
-QQ5 
1 55+0* 

-0.06 


j 32+0.94 
1 56+^ : ^ 

1 ' JO -0.36 


0.96 
1.08 


4 


84 +0 05 
U -°M>.04 


i 59+O.O4 
1-J -0.04 


1 56+ 029 


1.43 


2 20 +0 ' 27 
z ' zu -0.28 


1 03 +0 07 
■ UJ -0.08 


1 56 +0 07 

1 ' JO -0.06 


2 13+ - 32 

z - 1J -0.36 


1.15 


5 


0.84^ 


1 47+O.O4 
1 -^'-0.04 


1 70+O.I6 
1J -0.12 


1.35 


9 99+0.25 
z - z -0.25 


1 ns+o.06 


1 56^ m 

i- JD -0.12 


2 47+0.39 
z - H '-0.97 


1.32 


6 




1 49+0.04 
-0.04 


1 91+0.21 
1 " 41 -0.15 


1.55 


9 74+0.28 
z - -0.29 


1 10 +0 07 
1-1 -0.09 


| 47+O.O6 
1 - H '-0.05 


2 Ol^' 42 
Z - -0.42 


1.28 


7 


0.82+°$ 


i 4C+0.04 


1 7S+0-15 
'• zS -0.12 


1.05 


2.02+0-41 


097+0.13 


1 46+0- 05 


1 46+0- 46 

1 '^°-0.26 


1.00 


8-9 




| 44+0.03 


1 97+O.H 
'•^'-O.ll 


1.71 


2.00!"; 34 


1 01 +0 13 
1-ul -o.io 


| 44+O.O4 


1 41+0.75 
1 '^ 1 -0.17 


1.57 


10-12 


90+° 03 

u -' u -0.03 


| 47+O.O3 
1 -^'-0.03 


1 74+0-10 
1 - J ^-0.09 


0.96 




96 +0 08 

u -'°-0.09 


1 47+O.O3 
1 -+'-0.03 


1 41+O.I8 


0.93 


13-16 


95+0-03 

-0.03 


1 59+0.03 
1-J -0.03 


1 46+0- 16 


1.38 


9 95+0.21 
-0.20 


1 16+° 05 

1-1 -0.05 


1 58+ 0i)6 

• JO -0.06 


9 75+0.27 
Z - J -0.34 


1.01 


1998 April 23 rd 


1-5 
6-11 


96+° 03 
-0-Q3 
98 + <> 02 

u -'°-0.02 


1 58+ 003 

1 67 +ii 82 
1 - D '-0.02 


1.40+0.12 
1 7 f) +8:o9 

1,JU -0.07 


1.76 
2.45 


2 09+ - 32 

Z - U7 -Q.25 
, 99+0.21 


1 17+O.IO 

-QQ9 
1 09+008 

1 ' u ^-0.08 


1 59+0.07 

1 fi o+8:81 

1.68_ .0 3 


1.70+067 

1 43+^ : ^ 

-0.13 


1.50 
2.31 


12-19 


1 07+° 02 

^ -0.02 


1 69+° 03 


1 46+ 009 


1.79 


2 09+0.18 
Z - UZ -0.18 


1 20+ 006 


1 69 +0 03 


1 62 +0 19 


1.37 



a The Galactic hydrogen equivalent absorbing column is Nh = (1.61 ±0.1) X 10 20 cm 2 (Lockman & Savage 1995). The 
quoted errors are for 1 a for 3 parameters (i.e. A\ 2 = 3.53). The number of degrees of freedom is 57. 

b The quoted errors are for 1 a for 4 parameters (i.e. A\ 2 = 4.72). The number of degrees of freedom is 56. 
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